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ABSTRACT 

A debate has arisen regarding the importance of stationary versus eruptive 
mass loss for massive star evolution. The reason is that stellar winds have been 
found to be clumped, which results in the reduction of unclumped empirical mass- 
loss rates. Most stellar evolution models employ theoretical mass-loss rates which 
are already reduced by a moderate factor of ~2-3 compared to non-corrected em- 
pirical rates. A key question is whether these reduced rates are of the correct 
order of magnitude, or if they should be reduced even further, which would mean 
that the alternative of eruptive mass loss becomes necessary. Here we introduce 
the transition mass-loss rate M trans between O and Wolf-Rayet (WR) stars. Its 
novelty is that it is model independent. All that is required is postulating the 
spectroscopic transition point in a given data-set, and determining the stellar 
luminosity, which is far less model dependent than the mass-loss rate. The tran- 
sition mass-loss rate is subsequently used to calibrate stellar wind strength by its 
application to the Of/WNh stars in the Arches cluster. Good agreement is found 
with two alternative modelling/theoretical results, suggesting that the rates pro- 
vided by current theoretical models are of the right order of magnitude in the 
~5OM mass range. Our results do not confirm the specific need for eruptive 
mass loss as Luminous Blue Variables, and current stellar evolution modelling 
for Galactic massive stars seems sound. Mass loss through alternative mecha- 
nisms might still become necessary at lower masses, and/or metallicities, and the 
quantification of alternative mass loss is desirable. 

Subject headings: Stars: early-type - Stars: mass-loss - Stars: winds, outflows - 
Stars: evolution 



1. Introduction 



Mass loss via stellar winds is thought to play a dominant role in the evolution of massive 
O-type stars, because of the loss of mass, as winds "peel off" the star's outermost layers 
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(Conti 1976), as well as through the associated loss of angular momentum (e.g. Langer 
1998, Meynet & Maeder 2002). However, during the last decade, large uncertainty has been 
pointed out regarding our quantitative knowledge of the mass-loss rates of massive stars, as 
stellar winds have been revealed to be clumped, resulting in empirical rates that have been 
overestimated. 

Although it had been known for decades that O-type winds are clumped (Lupie & 
Nordsieck 1987, Eversberg et al. 1998), the severity did not appear to be fully recognized 
until Bouret et al. (2005) and Fullerton et al. (2006) claimed mass-loss reductions of factors 
~3-7 and ~20-130 respectively in comparison to unclumped Ha and radio mass-loss rates 
(e.g. Lamers & Leitherer 1993). The Ha diagnostics depends on the density squared, and 
are thus sensitive to clumping, whilst ultraviolet P Cygni lines such as Pv are insensitive to 
clumping as these depend linearly on the density. The above-mentioned Bouret et al. and 
Fullerton et al. analyses were based on models where the wind is divided into a portion of 
the wind containing all the material with a volume filling factor fy (the reciprocal of the 
clumping factor), whilst the remainder of the wind is assumed to be void. This pure micro- 
clumping approach is probably an oversimplification of the real situation, but it provides 
interesting insights into the potential mass-loss rate reductions. 

In reality, clumped winds are likely porous, with a range of clump sizes, masses, and 
optical depths. Macro-clumping and porosity have been investigated with respect to both 
the spectral analyses (e.g. Oskinova et al. 2007, Sunqvist et al. 2010, Surlan et al. 2012) 
as well as the radiative driving (Muijres et al. 2011). The upshot from these studies is 
that O star mass-loss rates may only be reduced by a moderate factor of ~3 (Repolust et 
al. 2004, Puis et al. 2008), which would bring their clumping properties in agreement with 
those of Wolf-Rayet (WR) winds, for which similar moderate clumping factors have been 
derived (Hamann & Koesterke 1998). The latter are based on the analysis of emission line 
wings due to electron scattering, which have the advantage that they do not depend on 
detailed ionization fractions and abundances of trace elements. These moderate clumping 
factors would imply that massive star evolution modelling is not affected, as current state- 
of-the-art rotating stellar models (e.g. Georgy et al. 2011; Brott et al. 2011) already employ 
moderately reduced rates via the theoretical relations of Vink et al. (2000). 

In light of the severe mass-loss reductions claimed e.g. by Fullerton et al. (2006), Smith 
& Owocki (2006) argued that the integrated mass loss from stationary stellar winds for very 
massive stars (VMS) above ~50M Q may be vastly insufficient to explain their role as the 
progenitors of WR stars and stripped-envelope Ibc supernovae. Instead, Smith & Owocki 
argued that the bulk of VMS mass loss is likely of an eruptive rather than a stationary 
nature. In particular, they highlighted the alternative option of eruptive mass loss during 
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the Luminous Blue Variable (LBV) phase. 

In view of the new porosity results, the arguments of Smith & Owocki (2006) however 
seem to have lost weight. Furthermore, quantitative estimates on the integrated amount of 
eruptive mass loss are hard to obtain as both the eruption frequency, and the amounts of 
mass lost per eruption span a wide range (of a factor 100) with LBV nebular mass estimates 
varying from ~O.1M in PCygni to ~10M Q in 77 Car, as discussed by Smith & Owocki 
(2006). Moreover, the energies required to produce such giant mass eruptions are very high 
(~ 10 50 erg), and their energy source is unknown. Soker (2004) discussed that the energy and 
angular momentum required for 7] Car great eruption cannot be explained within a single-star 
scenario. 

Whilst stationary winds in O and WR stars are ubiquitous, it is not at all clear if 
LBV-type objects like rj Car have encountered a special evolution (such as a merger) or if 
all massive stars go through eruptive mass-loss phases. On the other hand, for the most 
massive main-sequence WNh stars (Crowther et al. 2010, Bestenlehner et al. 2011) there is 
both theoretical and empirical evidence for strong Eddington parameter T-dependent (see 
definition Eq.[3]) mass loss (Grafener et al. 2011). For VMS the role of stationary mass loss 
has thus increased rather than decreased in recent years. 

In summary, the relevant roles of eruptive versus stationary mass loss seem rather un- 
certain and unsettled at the current time. There are ongoing debates as to whether wind 
clumping reduces the mass-loss rates by moderate factors of ~2-3, such that stellar evo- 
lution would not be affected, or by more severe factors of order ~10. In the latter case, 
line-driving would become negligible and alternatives such as eruptive mass loss would need 
to be considered. 

In order to address the relative role of wind versus eruptive mass loss, it would be 
beneficial to be able to calibrate either one of them. At the moment both stellar wind and 
eruptive mass loss could be inaccurate by factors of 10, and possibly even more. In this 
Letter, we attempt to alleviate this problem by presenting a methodology that involves a 
model-independent mass-loss indicator, the transition mass-loss rate M trans - located right 
at the transition from optically thin to optically thick stellar winds. Martins et al. (2008) 
found two mass-loss relations for VMS Arches cluster stars, one for the Of stars and one for 
the late-type WNh stars respectively. The fact that Wolf-Rayet stars with WNh spectral 
classification have optical depth larger than one has already been discussed in literature (e.g. 
Grafener & Hamann 2008), and one might thus expect to witness a transition from optically 
thin O-type winds to optically thick Wolf-Rayet winds. 

Vink et al. (2011) discovered a sudden change in the slope of the mass- loss versus T 



-4- 



relation at the transition from O-type (optically thin) to WR-type (optically thick) winds. 
Interestingly, this transition was found to occur for a wind efficiency parameter rj = M 
Voo/ (L/c) of order unity This key result from Monte Carlo modelling that the transition 
from O to WR-type mass loss coincides with rj ~ 1, can also be found analytically (Sect.|2]). 
And the result can be utilized to "calibrate" wind mass loss in an almost model independent 
manner (Sect. E}. 



2. The transition mass-loss rate 

Netzer & Elitzur (1993) and Lamers & Cassinelli (1999; hereafter LC99) give general 
momentum considerations for dust-driven winds (see LC99 pages 152-153) that can also be 
applied to line-driven winds. The integral form of the momentum equation contains four 
terms (Eq. 7.5 of LC99). Because hydrostatic equilibrium is a good approximation for the 
subsonic part of the wind, and the gas pressure gradient is small beyond the sonic point, 
LC99 argue that the second and third terms are negligible compared to the first and fourth, 
resulting in: 

°° „ 2 dv , f°° GM . . 2 
4 it y ' ■ ■ ' ^ \ \ , , \ ~~ . . — 



4nr 2 pv—dr + / — — (1 - T)pAnr z dr = 0. (1) 
dr L r z 



Employing the mass-continuity equation M = 4irr 2 pv, one obtains 
r°° . dv ■ r°° 

/ M—dr = Mv 00 =4:7iGM (T(r) - l)pdr (2) 
JR+ dr J rs 

where r s denotes the sonic radius and T(r) the Eddington factor with respect to the 
total flux-mean opacity ftp: 

Using the wind optical depth r = J r °° ttppdr, one obtains 

M Voo ~ r - l r = — r. (4) 

K C 1 

Where it is assumed that T is significantly larger than one, and the factor ^-j^- is thus 
close to unity (LC99's second assumption), resulting in 
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M VoD = -r. (5) 
c 

One can now derive a key condition for the wind efficiency number rj 

r / = T7 - = r = l. (6) 

The key point of our Letter is that we can employ the unique condition r\ — r = 1 right 
at the transition from optically thin O-star winds to optically-thick WR winds, and obtain 
a model-independent M. In other words, if we were to have an empirical data-set available 
that contains luminosity determinations for O and WR stars, we can obtain the transition 
mass-loss rate M trans simply by considering the transition luminosity Z-trans arid the terminal 
velocity representing the transition point from O to WR stars: 

M trans = ^ (7) 

VooC 

We note that this transition point can be obtained by purely spectroscopic means, 
independent of any assumptions regarding wind clumping. 



2.1. Testing the assumptions 

In the above analysis we made two assumptions that we wish to check with numerical 
tests involving sophisticated hydrodynamic wind models (from Grafener & Hamann 2008) 
and simpler /3-type velocity laws, commonly used in O /WR wind modelling. The results are 
compiled in Tab. 1. We first confirm LC99's first assumption through the comparison of rj, 
determined from M and i^, to the approximate rj' values as computed from the right-hand- 
side of integral in Eq.fU where rj'L/c = AitGM J r ° (r(r) — l)pdr. Evidently, the values of r] 
and r)' agree at the few percent level, and the first LC99 approximation is verified. 

Second, we investigate the assumption that the term E=i in Eq.HJis close to unity by 
numerical integration of r = up dr. We obtain a correction factor /, which we define by 



M Voc EE f-T. 

C 



With this definition Eq.|H] becomes 

Mv 



L/c < 9 > 
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Table 1: Wind models. If not stated otherwise, the adopted stellar parameters are T* = 35 kK 
and log(L/L ) = 6.0 (with M = 60 M ). 



p 


log(M) 

[Me] 

L yr 


Voo 

[la.] 


V 


rf 


1 


J 


WR22, M = 78.1 M , log(L/L ) 


= 6.3 






HYD 


-4.868 


979.5 


0.326 


0.315 


1 1 HQ 


U.Z 1 1 


1.0 


-4.868 


979.5 


0.326 


0.310 


1.160 


0.281 


1.0 


-4.607 


1785.0 


1.085 


1.062 


i.yoo 


U.ool 


Voo/V, 


2SC = 2.5 


-»• (r- l)/r ~ 2.5/3.5 = 


0.714 




0.5 


-4.905 


2289.6 


1.4 


1.385 


1.681 


0.833 


0.5 


-5.148 


2289.6 


0.8 


0.793 


0.956 


0.837 


0.5 


-5.750 


2289.6 


0.2 


0.200 


0.238 


0.840 


1.0 


-4.905 


2289.6 


1.4 


1.383 


1.938 


0.722 


1.0 


-5.148 


2289.6 


0.8 


0.793 


1.091 


0.733 


1.0 


-5.750 


2289.6 


0.2 


0.198 


0.280 


0.713 


1.5 


-4.905 


2289.6 


1.4 


1.386 


2.293 


0.610 


1.5 


-5.148 


2289.6 


0.8 


0.794 


1.299 


0.616 


1.5 


-5.750 


2289.6 


0.2 


0.199 


0.330 


0.605 


Voo/v, 


= 1-5 ■ 


-»■ (r- i)/r ~ 1.5/2.5 = 


0.600 




0.5 


-4.683 


1373.8 


1.4 


1.365 


Z.U4D 


U.D54 


0.5 


-4.926 


1373.8 


0.8 


0.780 


1.17a 


O.D to 


0.5 


-5.528 


1373.8 


0.2 


0.195 


0.301 


0.664 


1.0 


-4.683 


1373.8 


1.4 


1.362 


OK/1/1 
Z.D44 


u.ozy 


1.0 


-4.926 


1373.8 


0.8 


0.780 


1 

l.Ool 


u.ozz 


1.0 


-5.528 


1373.8 


0.2 


0.196 


U.ooo 


U.OlD 


1.5 


-4.683 


1373.8 


1.4 


1.364 


3.231 


0.433 


1.5 


-4.926 


1373.8 


0.8 


0.779 


1 Q71 

i.y 1 1 


U.4UD 


1.5 


-5.528 


1373.8 


0.2 


0.196 


0.490 


0.408 


Voo/v, 


JSC = 2.5, 


= 70 kK 








0.5 


-5.449 


4579.2 


0.8 


0.800 


0.955 


0.837 


1.0 


-5.449 


4579.2 


0.8 


0.800 


1.140 


0.701 


1.5 


-5.449 


4579.2 


0.8 


0.800 


1.442 


0.555 


Voo/v, 


3SC = 2.5, = 17.5 kK 








0.5 


-4.847 


1144.8 


0.8 


0.784 


0.946 


0.846 


1.0 


-4.847 


1144.8 


0.8 


0.784 


1.033 


0.774 


1.5 


-4.847 


1144.8 


0.8 


0.783 


1.203 


0.665 


Voo/v, 


3SC = 2.5, log(L/L ) = 6.3 (with M = 120 M ) 


0.5 


-4.923 


2724.4 


0.8 


0.794 


0.959 


0.833 


1.0 


-4.923 


2724.4 


0.8 


0.794 


1.121 


0.714 


1.5 


-4.923 


2724.4 


0.8 


0.795 


1.366 


0.585 


Voo/v, 


3SC = 2.5, log(L/L ; 


) = 5.7 (with M = 40 M ) 


0.5 


-5.438 


2237.6 


0.8 


0.794 


0.953 


0.839 


1.0 


-5.438 


2237.6 


0.8 


0.794 


1.085 


0.737 


1.5 


-5.438 


2237.6 


0.8 


0.794 


1.276 


0.627 
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V ~ = fi'rad - -— —, (11) 



To compute the integral numerically, we need to obtain the density p(r), and the flux- 
mean opacity «p(r) in the stellar wind (T(r) follows from Eq. fl3]). The hydrodynamic wind 
models of Grafener & Hamann (2008) have these quantities directly available. We have 
performed a direct computation for the first model (HYD) in Tab. 1 for the Galactic WNh 
star WR22 (Grafener & Hamann 2008). For this model we obtain / = 0.272. This value is 
lower than, but of the order of, unity. The terminal wind speed in this model is significantly 
lower than the observed value for WR22 (980km/s vs. 1785 km/s). Consequently, our 
derived / is likely on the low side. We expect T = g Ta d/g to be connected to the ratio 
(foo + v esc )/v esc = foo/fesc + 1, and / ~ 

Here we follow a model-independent approach, adopting /3-type velocity laws. The mean 
opacity Kf then follows from the resulting radiative acceleration g ra d 

£rad(r) =<t f (r)^. (10) 

g ra d follows from the prescribed density p(r) and velocity structures v(r) via the equation of 
motion 

dv _ 1 dp GM 

'~T~ S'rad 3 n 

dr p dr r z 

where we assume a grey temperature structure to compute the gas pressure p. We note 
that these results are completely independent of any assumptions regarding wind porosity, 
or the chemical composition of the wind material. The only assumption that goes into 
these considerations is that the winds are radiatively driven. The resulting mean opacity Kp 
consequently captures all physical effects that could potentially affect the radiative driving. 

The obtained values for the correction factor / are summarized in Tab. 1. The first 
three models in Tab. 1 represent a consistency test with the hydrodynamic model for WR22. 
Using a beta law with = 1, and the same Voo, we obtain almost exactly the same / as for 
the hydrodynamic model, justifying our /3-law approach, which we employ in the following. 
Now employing the observed - and therefore likely close to correct - value of (and a 
correspondingly increased M, we obtain / = 0.55. 

To get a handle on the overall behaviour of this factor /, we computed a series of wind 
models for a range of stellar parameters 17.5 < T*/kK < 70 and 5.7 < log(L/L ) < 6.3, with 
wind efficiencies around the transition region (0.2 < rj < 1.4). Remarkably, the resulting 
values of / depend only on the adopted values of foo/^esc and (3. For foo/fesc = 2.5, we 
obtain / ~ 0.8, 0.7, 0.6 respectively for (3 = 0.5, 1.0, 1.5, where the last value is probably 
most appropriate (Vink et al. 2011). Overall, we derive values of / in the range 0.4-0.8, 
with a mean value of 0.6. We note that the error on this number / is within the uncertainty 
of the luminosity determinations described in the next section. 
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For transition objects with r ~ 1 we thus expect that / ~ 0.6, i.e. the transition between 
O and WR spectral types should occur at mass-loss rates of 

M = ~ 0.6M trans . (12) 

VooC 

The fact that the correction factor is within a factor of two of our idealized approach (Eq. [7]) 
is highly encouraging. We stress that this number is independent of any potential model 
deficiencies, as we have used the observed values of Vqq in this analysis. 

3. The transition mass loss rate in the Arches cluster 

Martins et al. (2008) analyzed 28 VMS in the Arches cluster, with equal numbers of 
O-type supergiants and nitrogen-rich Wolf-Rayet (WNh) stars (sometimes called "O stars on 
steroids"). For the O-type supergiants, we expect the winds to be optically thin, whilst the 
WNh stars should have optically thick winds. Here we postulate that the 04-6If+ represent 
the transition point where the optical depth crosses unity. 

In Tab. 2, we compiled a subset of 20 stars, skipping those objects with a He-enriched 
surface composition. The objects are sorted with respect to their spectral subtypes, and 
within each subtype bin with respect to their luminosity. Together with the basic stellar 
and wind parameters derived by Martins et al. (2008), we list the mass-loss rate for which 
M n= i if the stars would have a wind efficiency of exactly 1. The values listed in Tab. 2 show 
that there is a transition between O and WR spectral types. The spectroscopic transition 
for spectral subtypes 04-6If+ occurs at log(L) = 6.05 and \og(M v= i / MqJt' 1 ) = —4.95. 

This is the resulting transition mass-loss rate for the Arches cluster stars. Its determined 
value does not depend on model uncertainties involving issues such as wind clumping. The 
only remaining uncertainties are due to uncertainties in the terminal velocity and the stellar 
luminosity L. The latter results from errors in the distances and reddening parameters, as 
well as the determination of effective temperatures from non-LTE model atmospheres. If 
the derived value for M trans is compared with empirically determined mass-loss rates, the 
uncertainties in distances and reddening parameters nearly cancel, as for empirical mass-loss 
rates based on recombination line analyses M oc L 3 / 4 , while Mtrans oc L. 

To estimate uncertainties in the effective temperature scale for O stars we can use 
historical values from the last four decades (e.g. Panagia et al. 1973, Martins et al. 2005) as 
an indicator for potential systematic errors in the inclusion/neglect of certain micro-physics 
(line blanketing, wind effects, etc.), the best error estimate is ~10% in effective temperature, 
leading to potential errors in the luminosity of at most ~40%. This is several factors smaller 
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Table 2: The transition from O to WR stars for the most massive stars in the Arches Cluster. 



Star 


subtype 


log(L) 
[Le\ 


log(M) 
[^] 


Voo 
[knr] 


log(M„=i) 


F9 


WN8-9 


6.35 


-4.78 


1800 


-4.60 


Fl 


WN8-9 


6.30 


-4.70 


1400 


-4.54 


F14 


WN8-9 


6.00 


-5.00 


1400 


-4.84 


Bl 


WN8-9 


5.95 


-5.00 


1600 


-4.95 


F16 


WN8-9 


5.90 


-5.11 


1400 


-4.94 


F15 


04-6If+ 


6.15 


-5.10 


2400 


-4.92 


F10 


04-6If+ 


5.95 


-5.30 


1600 


-4 95 M 

^ ■■~" J lvl trans 


F18 


04-61 


6.05 


-5.35 


2150 


-4.98 


F21 


04-61 


5.95 


-5.49 


2200 


-5.09 


F28 


04-61 


5.95 


-5.70 


2750 


-5.18 


F20 


04-61 


5.90 


-5.42 


2850 


-5.25 


F26 


04-61 


5.85 


-5.73 


2600 


-5.26 


F32 


04-61 


5.85 


-5.90 


2400 


-5.22 


F33 


04-61 


5.85 


-5.73 


2600 


-5.26 


F22 


04-61 


5.80 


-5.70 


1900 


-5.17 


F23 


04-61 


5.80 


-5.65 


1900 


-5.17 


F29 


04-61 


5.75 


-5.60 


2900 


-5.41 


F34 


04-61 


5.75 


-5.77 


1750 


-5.19 


F40 


04-61 


5.75 


-5.75 


2450 


-5.33 


F35 


04-61 


5.70 


-5.76 


2150 


-5.33 



Notes - Designations, subtypes, luminosities (L), mass- loss rates (M), and terminal wind 
velocities (i>oo) according to Martins et al. (2008). The 6th column indicates the 
mass-loss rate where r]—l. For the Arches cluster, we obtain log(M trans ) ~ —4.95. 
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than the order-of-magnitude uncertainties in mass-loss rates due to clumping and porosity. 
In other words our simple equation (EqJTJ) is of tremendous value in calibrating stellar wind 
mass loss, and assessing its role in the mass loss during the evolution of massive stars (Sect.HJ). 

How does our determination of the transition mass-loss rate compare to other methods? 
Let us first compare our transition mass-loss value to the mass-loss rates of Martins et 
al. Martins et al. use the non-LTE cmfgen code by Hillier & Miller (1998), employing 
a micro-clumping approach with a volume filling factor / = 0.1 for their K-band analysis. 
Their values are in good agreement with our transition mass-loss rates for the objects at 
the boundary between O and WR (see their Table 2). This is unlikely to be a coincidence. 
According to our findings in Sect. 12. 1( we expect mass-loss rates of the order of M ~ 0.6 x 
^trans for the transition objects, i.e. log(M/M yr _1 ) = —5.2. 

We also compare the transition mass-loss rate to the oft-used theoretical mass-loss 
relation of Vink et al. (2000), for which we find logMymk = — 5.14M yr _1 for an assumed 
stellar mass M = 60 M . This number is within 0.2 dex from the transition mass-loss rate 
logM trans = — 4.95M yr _1 . The comparison is hardly compromised as a result of the Vink 
et al. dependence on stellar mass, as for masses in the range 40-80 M , the Vink et al. 
mass-loss rate varies by at most 0.04 dex. 

In summary, we have three independent mass-loss rate determinations that agree within 
a factor of two. This means that our concept of the transition mass-loss rate has indeed been 
able to test the accuracy of current mass-loss estimates by stellar winds. 

4. Discussion 

Now that we have calibrated stellar wind mass loss in the high mass and luminosity 
regime, we assess the role of stellar wind mass loss for massive star evolution. For a 60 M 
star, the main sequence lifetime is ~3 Myrs (e.g. Weidner & Vink 2010 and references 
therein). With a stationary mass-loss rate of 10 -5 M yr _1 as derived for the transition 
mass-loss rate in Sect. [3], this means such an object will lose ~30 M , i.e. half its initial mass, 
already on the main sequence during core hydrogen burning. We have not yet addressed the 
new concept of Y dependent mass loss, nor any additional stellar wind mass loss during the 
subsequent core-helium burning WR phase. In other words, solar-metallicity 60 M stars 
are expected to lose the bulk of their initial masses through stellar winds, leaving very little 
(if any) space for additional (e.g. eruptive) mass loss. 

It is plausible that the strong-winded VMS remain on the blue side of the HR diagram, 
without ever entering an eruptive Luminous Blue Variable (LBV) or Red supergiant (RSG) 
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phase. Current wisdom thus suggests that solar metallicity VMS likely "evaporate" primarily 
as the result of stationary wind mass loss, without the necessity of additional eruptive mass 
loss. However, finding out if eruptive mass loss might play an additional role remains an 
interesting exercise, especially for the lower initial mass and sub-solar metallicity ranges, 
as their story might be expected to be different. Moreover, we know 77 Car analogs and 
supernova impostors exist in external galaxies (e.g. Van Dyck et al. 2005, Pastorello et al. 
2010, Kochanek et al. 2012). 

Contrary to the most massive stars, stars below ~4OM likely evolve into the RSG, 
yellow super/hypergaint, and/or LBV regimes of the stellar HR diagram (see e.g. Vink 
2009). We note that the Vink et al. (2000) main-sequence mass- loss rates currently in use 
in stellar models (e.g. Brott et al. 2011) for lower mass "normal" 20-60 M & O stars are 
also already reduced by a factor 2-3 in comparison to previous unclumped empirical rates. 
There is currently no particular reason to assume they are still overestimated^. The results 
presented here certainly boost confidence in the mass-loss rates currently in use, although 
they remain uncalibrated for the lower mass regime. One should also realize when working 
down the mass range, starting from our 6OM calibrator star, that the mass-loss rates drop 
significantly below 10~ 5 M yr _1 . Its effects on stellar evolution remain significant due to the 
longer evolutionary timescales for lower mass objects and the fact that it is the multiplication 
of the mass- loss rate times the duration that is relevant. This is especially relevant for angular 
momentum evolution, possibly down to stellar masses as low as 10-15 M (Vink et al. 2010). 
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